Comparison of stochastic and deterministic models for gambiense sleeping sickness at different spatial scales: A health area analysis in the DRC.

The intensification of intervention activities against the fatal vector-borne disease gambiense human African trypanosomiasis (gHAT, sleeping sickness) in the last two decades has led to a large decline in the number of annually reported cases. However, while we move closer to achieving the ambitious target of elimination of transmission (EoT) to humans, pockets of infection remain, and it becomes increasingly important to quantitatively assess if different regions are on track for elimination, and where intervention efforts should be focused. We present a previously developed stochastic mathematical model for gHAT in the Democratic Republic of Congo (DRC) and show that this same formulation is able to capture the dynamics of gHAT observed at the health area level (approximately 10,000 people). This analysis was the first time any stochastic gHAT model has been fitted directly to case data and allows us to better quantify the uncertainty in our results. The analysis focuses on utilising a particle filter Markov chain Monte Carlo (MCMC) methodology to fit the model to the data from 16 health areas of Mosango health zone in Kwilu province as a case study. The spatial heterogeneity in cases is reflected in modelling results, where we predict that under the current intervention strategies, the health area of Kinzamba II, which has approximately one third of the health zone's cases, will have the latest expected year for EoT. We find that fitting the analogous deterministic version of the gHAT model using MCMC has substantially faster computation times than fitting the stochastic model using pMCMC, but produces virtually indistinguishable posterior parameterisation. This suggests that expanding health area fitting, to cover more of the DRC, should be done with deterministic fits for efficiency, but with stochastic projections used to capture both the parameter and stochastic variation in case reporting and elimination year estimations.

Trends in outcomes following COVID-19 symptom onset in Milan: a cohort study.

BACKGROUND: For people with symptomatic COVID-19, the relative risks of hospital admission, death without hospital admission and recovery without admission, and the times to those events, are not well understood. We describe how these quantities varied with individual characteristics, and through the first wave of the pandemic, in Milan, Italy. METHODS: A cohort study of 27 598 people with known COVID-19 symptom onset date in Milan, Italy, testing positive between February and June 2020 and followed up until 17 July 2020. The probabilities of different events, and the times to events, were estimated using a mixture multistate model. RESULTS: The risk of death without hospital admission was higher in March and April (for non-care home residents, 6%-8% compared with 2%-3% in other months) and substantially higher for care home residents (22%-29% in March). For all groups, the probabilities of hospitalisation decreased from February to June. The probabilities of hospitalisation also increased with age, and were higher for men, substantially lower for healthcare workers and care home residents, and higher for people with comorbidities. Times to hospitalisation and confirmed recovery also decreased throughout the first wave. Combining these results with our previously developed model for events following hospitalisation, the overall symptomatic case fatality risk was 15.8% (15.4%-16.2%). CONCLUSIONS: The highest risks of death before hospital admission coincided with periods of severe burden on the healthcare system in Lombardy. Outcomes for care home residents were particularly poor. Outcomes improved as the first wave waned, community healthcare resources were reinforced and testing became more widely available.

Risk factors associated with severe hospital burden of COVID-19 disease in Regione Lombardia: a cohort study.

BACKGROUND: Understanding the risk factors associated with hospital burden of COVID-19 is crucial for healthcare planning for any future waves of infection. METHODS: An observational cohort study is performed, using data on all PCR-confirmed cases of COVID-19 in Regione Lombardia, Italy, during the first wave of infection from February-June 2020. A multi-state modelling approach is used to simultaneously estimate risks of progression through hospital to final outcomes of either death or discharge, by pathway (via critical care or not) and the times to final events (lengths of stay). Logistic and time-to-event regressions are used to quantify the association of patient and population characteristics with the risks of hospital outcomes and lengths of stay respectively. RESULTS: Risks of severe outcomes such as ICU admission and mortality have decreased with month of admission (for example, the odds ratio of ICU admission in June vs March is 0.247 [0.120-0.508]) and increased with age (odds ratio of ICU admission in 45-65 vs 65 + age group is 0.286 [0.201-0.406]). Care home residents aged 65 + are associated with increased risk of hospital mortality and decreased risk of ICU admission. Being a healthcare worker appears to have a protective association with mortality risk (odds ratio of ICU mortality is 0.254 [0.143-0.453] relative to non-healthcare workers) and length of stay. Lengths of stay decrease with month of admission for survivors, but do not appear to vary with month for non-survivors. CONCLUSIONS: Improvements in clinical knowledge, treatment, patient and hospital management and public health surveillance, together with the waning of the first wave after the first lockdown, are hypothesised to have contributed to the reduced risks and lengths of stay over time.

Decreasing hospital burden of COVID-19 during the first wave in Regione Lombardia: an emergency measures context.

BACKGROUND: The aim of this study is to quantify the hospital burden of COVID-19 during the first wave and how it changed over calendar time; to interpret the results in light of the emergency measures introduced to manage the strain on secondary healthcare. METHODS: This is a cohort study of hospitalised confirmed cases of COVID-19 admitted from February-June 2020 and followed up till 17th July 2020, analysed using a mixture multi-state model. All hospital patients with confirmed COVID-19 disease in Regione Lombardia were involved, admitted from February-June 2020, with non-missing hospital of admission and non-missing admission date. RESULTS: The cohort consists of 40,550 patients hospitalised during the first wave. These patients had a median age of 69 (interquartile range 56-80) and were more likely to be men (60%) than women (40%). The hospital-fatality risk, averaged over all pathways through hospital, was 27.5% (95% CI 27.1-28.0%); and steadily decreased from 34.6% (32.5-36.6%) in February to 7.6% (6.3-10.6%) in June. Among surviving patients, median length of stay in hospital was 11.8 (11.6-12.3) days, compared to 8.1 (7.8-8.5) days in non-survivors. Averaged over final outcomes, median length of stay in hospital decreased from 21.4 (20.5-22.8) days in February to 5.2 (4.7-5.8) days in June. CONCLUSIONS: The hospital burden, in terms of both risks of poor outcomes and lengths of stay in hospital, has been demonstrated to have decreased over the months of the first wave, perhaps reflecting improved treatment and management of COVID-19 cases, as well as reduced burden as the first wave waned. The quantified burden allows for planning of hospital beds needed for current and future waves of SARS-CoV-2 i.

Forecasting the 2017/2018 seasonal influenza epidemic in England using multiple dynamic transmission models: a case study.

BACKGROUND: Since the 2009 A/H1N1 pandemic, Public Health England have developed a suite of real-time statistical models utilising enhanced pandemic surveillance data to nowcast and forecast a future pandemic. Their ability to track seasonal influenza and predict heightened winter healthcare burden in the light of high activity in Australia in 2017 was untested. METHODS: Four transmission models were used in forecasting the 2017/2018 seasonal influenza epidemic in England: a stratified primary care model using daily, region-specific, counts and virological swab positivity of influenza-like illness consultations in general practice (GP); a strain-specific (SS) model using weekly, national GP ILI and virological data; an intensive care model (ICU) using reports of ICU influenza admissions; and a synthesis model that included all data sources. For the first 12 weeks of 2018, each model was applied to the latest data to provide estimates of epidemic parameters and short-term influenza forecasts. The added value of pre-season population susceptibility data was explored. RESULTS: The combined results provided valuable nowcasts of the state of the epidemic. Short-term predictions of burden on primary and secondary health services were initially highly variable before reaching consensus beyond the observed peaks in activity between weeks 3-4 of 2018. Estimates for R0 were consistent over time for three of the four models until week 12 of 2018, and there was consistency in the estimation of R0 across the SPC and SS models, and in the ICU attack rates estimated by the ICU and the synthesis model. Estimation and predictions varied according to the assumed levels of pre-season immunity. CONCLUSIONS: This exercise successfully applied a range of pandemic models to seasonal influenza. Forecasting early in the season remains challenging but represents a crucially important activity to inform planning. Improved knowledge of pre-existing levels of immunity would be valuable.

Exploiting routinely collected severe case data to monitor and predict influenza outbreaks.

BACKGROUND: Influenza remains a significant burden on health systems. Effective responses rely on the timely understanding of the magnitude and the evolution of an outbreak. For monitoring purposes, data on severe cases of influenza in England are reported weekly to Public Health England. These data are both readily available and have the potential to provide valuable information to estimate and predict the key transmission features of seasonal and pandemic influenza. METHODS: We propose an epidemic model that links the underlying unobserved influenza transmission process to data on severe influenza cases. Within a Bayesian framework, we infer retrospectively the parameters of the epidemic model for each seasonal outbreak from 2012 to 2015, including: the effective reproduction number; the initial susceptibility; the probability of admission to intensive care given infection; and the effect of school closure on transmission. The model is also implemented in real time to assess whether early forecasting of the number of admissions to intensive care is possible. RESULTS: Our model of admissions data allows reconstruction of the underlying transmission dynamics revealing: increased transmission during the season 2013/14 and a noticeable effect of the Christmas school holiday on disease spread during seasons 2012/13 and 2014/15. When information on the initial immunity of the population is available, forecasts of the number of admissions to intensive care can be substantially improved. CONCLUSION: Readily available severe case data can be effectively used to estimate epidemiological characteristics and to predict the evolution of an epidemic, crucially allowing real-time monitoring of the transmission and severity of the outbreak.